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ABSTRACT 

We made a narr owband NB973 (bandwidth of 200A centered at 9755A) imaging of the Subaru Deep 
Field (SDF) and found two z — 7 Lya emitter (LAE) candidates down to NB973 = 24.9. Carrying out 
deep follow-up spectroscopy, we identified one of them as a real z — 6.96 LAE. This has established a 
new redshift record, showing that galaxy formation was in progress just 750 Myr after the Big Bang. 
Meanwhile, the Lya line luminosity function of LAEs is known to decline from z — 5.7 to 6.6 in the 
SDF. L* at z = 6.6 is 40-60% of that at z = 5.7. We also confirm that the number density of z = 7 
LAE is even only 17% of the density at z = 6.6 comparing the latest SDF LAE samples. This series of 
significant decreases in LAE density with increasing redshift can be the result of galaxy evolution during 
these epochs. However, using the UV continuum luminosity functions of LAEs, those of Lyman break 
galaxies and a LAE evolution model based on hierarchical clustering, we find that galaxy evolution alone 
cannot explain all the decrease in density. This extra density deficit might reflect the attenuation of the 
Lya photons from LAEs by the neutral hydrogen possibly left at the last stage of the cosmic reionization 
at z ~ 6-7. 

Subject headings: cosmology: observations — early universe — galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

Investigating high redshift galaxies as well as other dis- 
tant objects in the early Universe, especially within the 
first 1 Gyr after the Big Bang, has been the key to un- 
derstanding how galaxies have formed and evolved, probe 
their star formation histories, and constrain the epoch 
of cosmic reionization. The latest measurements of the 
polarization of the cosmic microwave background (CMB) 
by Wilkinson Microwave Anisotropy Probe ( WMAP) con- 
strained the optical depth to electron scattering dur- 
ing reionization and suggests that the average redshift 
of reionization was z = 10.9t2:3 ('Spcrgel et all l2007t 
iPage et all I2007D. Also, Gunn-Peterson (CP) troughs 
(iGunn &: PetersonI [19651) in z ~ 6 quasar spectra im- 
ply reionization ended at z ~ 6 with an estimated frac- 
tion of intergalactic medium (IGM ) neutral hydrogen, 
xgj'6.2 ^ 0.01-0.04 (|Fan et al.|[200l) . Moreover, a spec- 
tral modeling analysis of a z ~ 6.3 gamma ray burst 
(GRB) shows that the Universe seems to have been largely 
reionized at z ~ 6.3 with x^^-^ = and the upper 
limit of x^^'^ < 0.17-0.6, which suggests only some rea- 
sonable amount of neutral gas in the GRB host galaxy 
(jTotani et al.ll2006t) . 

Another probe of reionization is Lya emitters (LAEs), 
young galaxies in the distant universe showing in their 



spectra redshifted Lya emission from their interstellar gas 
illuminated by massive stars. The observed Lya line lu- 
minosity function (Lya LP) is expected to decline beyond 
the redshift z ~ 6 where reionization is thought to have 
completed as the increasing fraction of IGM neutral hy- 
drogen a bsorbs or scatters the Lya photons from young 
galaxies (iHainian fc Spaansl 119991 : iRhoads fc Malhotra I 
1200 It iHu et al.l l2002f). Nevertheless, recent LAE surveys 
show that Lya LP seems not to change fr om z = 3 to 
5.7 dOuchi et al. 2003; Aiiki et al. 2003: T r an et al.l 12004 
Ivan Breukelen et a l. 2005; Ouchi et al.ll2007l) . For the ear- 
lier epoch, Malhotra fc Rhoadsl ( 20041 ) suggest that Lya 
LP does not evolve between z = 5.7 and 6.6. This might 
be because their sample could be somewhat biased since 
it consists of several LAE subsamples taken from various 
surveys with different kinds of factors such as selection 
criteria, analysis methods, sky areas, survey volumes, and 
depths in order to compile as large a sample as possible. 

On the other ha nd, the Subaru Deep Field 
(|Kashikawa et al.ll2004 SDF) surveys have tried to keep 
all these factors as consistent as possible among differ- 
ent redshifts, surveyed exceptionally large volume and 
made large amount of LAE samples aX z — 4.8, 5.7 and 
6.6. Their latest survey has for the first time confirmed 
that the Lya LP declines as i^^g.e ^ -^z=5.7>< (0-4-0.6) 
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from z = 5.7 to 6.6 even after correcting cosmic variance 
(jKashikawa et al.ll2006H ). From this decline of the LF, 
they estimated the upper hmit of the neutral fraction at 
z = 6.6 to be < xgf < 0.45. If the nenutral IGM re- 
mains at ~ 50% level at z = 6.6, this constraint supports 
late reionization and contradicts the WMAP result. Also, 
the decline of Lya LF at z = 5.7-6.6 can also be ascribed 
to the evolution of LAE population itself. 

Meanwhile, the ionized fraction Xi < 1 and the mor- 
phology of Hii regions during patchy reionization would 
modulate the observed distribution of LAEs and enhance 



the observed cluster i ng of them (iFurlanetto et al.l l2006t 
iMcQuinn et al.ll2007D . iMcQuinn et al.l (|2007f) investigated 



the angular correlation function of the SDF photomet- 
ric sam ple of z = 6.6 LAEs obtained by iKashikawa et "al] 
(|2006bD and suggest that the Universe is fully ionized at 
z = 6.6 with th e mean volume ionized fraction of af^ ~ 1. 
IMcQuinn et al.l (|2007D also pointed out the difficulty in 
distinguishing the effect of evolution of LAE population 
on Lya LF from that of reionization. 

LFs of high-z galaxies also tell us about galaxy evolution 
itself in terms of how many galaxies existed at each lumi- 
nosity and epoch in the history of the Universe and how 
it has changed with cosmic time. To obtain this kind of 
information at high redshifts, LFs of Lyman break galax- 
ies (LBGs) and LAEs have been mainly observed. Ultra- 
violet continuum luminosity functions (UVLFs) of LBGs 
have been investigated from z ^ 3 to z ^ 7 and found 
to decline as redshift increases (iLehnert fc Bre mer "2003!: 



Ouchi ct al.1 |2004 i Bouwc ns et al.l 120061 : lYoshida ct al 



20061: iBouwens k Illingworth[|2006D . Since the UV contin 



uum redder than 1216A is not attenuated by neutral IGM 
hydrogen and if dust extinction is precisely corrected, the 
decline of UVLF reflects the evolution of galaxies. 

One recently observed example of this is a large decline 
of UVLF of dropout galaxies at 6 < z < 7-8, which is con- 
sidered to be a clear sign of galaxy ev olution over these 
redshifts (jBouwens fc IllingworthI [20061 ) . They conclude 
that very luminous galaxies are quite rare at z = 7-8. On 
the other hand, the U VLF of LAEs was confirmed not 
to evolve ai z ^ 3-5 (jOuchi et al.l |2003[) . In addition, 
studying LAEs in an even wider sky region, ~ 1.0 deg^ 
of the Suh& TVLl XMM-Newton Deep Survey (SXDS) field, 
lOuchi et all (|2007l ) found that UVLF of LAEs increases 
from z~3 — 4to5.7 while Lya LF of them remains un- 
changed over these redshifts, suggesting that the fraction 
of UV-bright LAEs increases at z = 5.7. Furthermore, no 
evolution of LAE UVLF from z = 5.7 to 6.6 was observed 
while Lya LF of LAE s evolves between these epochs in 
the latest SDF survey (jKashikawa et al.ll2006"H) . This im- 
plies that LAEs themselves do not significantly evolve from 
z = 6.6 to 5.7 and the decline of the Lya LF might reflect 
the effect of reionization. However, we do not know if this 
trend of the Lya LF and UVLF of LAEs continues from 
even earlier epochs. In other words, it is not clear whether 
LAE population evolves from z > 6.6 as LBG does and 
the neutral IGM fraction increases to suppress the Lya 
LF more severely beyond z — 6.6. Moreover, the existence 
of the galaxies at z > 6.6 has not been confirmed by spec- 
troscopy yet though several photometric candidates have 
been found. These questions can be addressed by observ- 
ing LAEs and their LFs at z > 6.6. Investigating their 



change over longer cosmic time interval, we can constrain 
the galaxy evolution and reionization more tightly. 

One possible method of detecting z > 6.6 LAEs is a 
narrowband filter imaging in infrared region. However, be- 
yond the redward limit of CCD sensitivity, the large format 
mosaicing advantages of infrared arrays are not yet avail- 
able and observations of high redshift LAEs is limited to 
a small survey volume. Though recent infrared detectors 
have achieved extremely high sensitivities, surveys with 
them cannot avoid large uncertianty due to cosmic vari- 
ance. Therefore, we carried out a narrowband survey of 
z = 7 LAEs using the final window of OH-airglow at the 
very edge of the optical regime still accessible with CCDs 
of Subaru Prime Focus Camera (jMivazaki et al.l l2002l 
Suprime-Cam) having a superb wide field of view, 34' x 27'. 
We chose the wavelength region 9655-9855A open to the 
highest redshift optical narrowband survey. Although this 
might not be considered quite an adequate window since 
there are several OH lines in the region, the estimated frac- 
tion of the sky counts coming from OH lines in the win- 
dow is not prohibitively large (only ~ 4.3 photons s~^ A"-'^ 
arcsec^^ m^^) and we actually succeeded in making a nar- 
rowband filter named NB973 covering this wavelength re- 
gion. This range corresponds to a redshift of 6.9 < z < 7.1 
for LAEs. 

To discover such extremely high redshift LAEs and make 
a sample of them as consistent a s possible with those o f 
z = 5.7 and 6.6 LAEs obtaine d bv lShimasaku et all (|2006f) 
and ' Kashikawa et al.l (j2006b[ ). we targeted the same field, 
SDF using 8.2-m Subaru Telescope/Suprime-Cam. Our 
brief and preliminary result has been recently reported in 
llveet all (120061 In this survey, we successfully confirmed 
a z = 6.96 LAE spectroscopically and observed that the 
number density of z = 7 LAE further declines from z = 6.6 
by a factor of 0.18-0.36, suggesting that the neutral hy- 
drogen might increase between these epochs. However, we 
do not know whether there had been any possible evolu- 
tion of LAE population itself from z = 7 to 6.6, and the 
density dificit might come from such a galaxy evolution. 

In this paper, we present the methods and results of our 
photometric and spectroscopic surveys fo r z = 7 LAEs, 
which were not fully covered in live et al.l (|2006l ). and try 
to draw out as much useful information as possible about 
the epoch of reionization and the LAE galaxy evolution 
from our results combined with the most recent high red- 
shift galaxy surveys and a LAE evolution model based on 
hierarchical clustering. We first describe the properties 
of our narrowband filter NB973 and imaging observation 
in §2. Then, selection criteria of z = 7 LAE candidates 
based on narrow and broadband images are derived in de- 
tail and their photometric properties are analyzed in §3. In 
§4, we explain the results of our follow-up spectroscopy of 
the selected candidates and their spectroscopic properties. 
In §5, we compare Lya and UV LFs of z = 7 LAEs with 
those of z = 5.7 and 6.6 LAEs derived fro m the latest sam- 
ples obtained by the SDF LA E surveys (jShimasaku et al.l 
120061: IKashikawa et al.l l2006b[ ) and discuss what implica- 
tions the result gives for cosmic reionization and galaxy 
evolution. Also, any possibilities of LAE galaxy evolution 
at z = 5.7-7 are inspected by observational and theor- 
eical approaches. In the last section, we conclude and 
summarize our results. Throughout we adapt a concor- 
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dance cosmology with {ilm, ^A, h) = (0.3,0.7,0.7), and 
AB magnitudes with 2" diameter aperture unless other- 
wise specified. 

2. IMAGING OBSERVATION 

We developed a narrowband filter, NB973, designed to 
cover the last optical window of OH-airglow centered at 
9755A with AAfwhm ~ 200A corres ponding to Lya emis- 
sion at 6.9 < z < 7.1 (|Ive et al.ll2006t ). The design and fab- 
rication of such a narrowband filter was not a simple issue 
for Suprime-Cam that uses a fast converging F/1.83 beam 
whose incident angle varies with a position in the field of 
view. Mixture of such light with different incident angles 
severely degrades the resultant transmission characteris- 
tics of the narrowband filter from our target design. To 
obtain the desired performance complying with the filter 
specification for our scientific requirement, we employed a 
combination of the following three filters glued together: 
a color cut glass filter RG780 with anti-reflection coating, 
a narrow bandpass multi-layer coating filter, and another 
multi-layer coating filter for red leak prevention. One year 
before starting the manufacturing of the NB973 filter for 
Suprime-Cam, we made another filter NB980 (bandwidth 
of - lOOA centered at - 9800A) for use in the paral- 
lel beam section of the Faint Ob ject Camera And Spec- 
trograph (jKashikawa et al.l |2003 . FOCAS) on Subaru to 
demonstrate the feasibility of narrowband imaging at this 
last OH window. During this prefabrication of NB980, 
manufacturing errors in controlling the thickness of thin 
film layers were evaluated. The multilayer thin film coat- 
ing design for NB973 was then optimized so that the re- 
sulting transmitting properties are relatively robust to in- 
evitable manufacturing errors to control the thickness of 
each thin layers. The measured transmission curve of the 
final NB973 filter actually used in the present survey as 
well as other filters used for color selection of z = 7 LAE 
candidates are shown in Figure [1] 

Our target sky region is SDF (jKashikawa et all |2004 
13^24'^21.M, -27°29'23"(J2000), ~876 arcmin^), a blank 
field in which z = 5.7 and 6.6 LAE surveys ha d been 
also ca rried out (jShimasaku et a l. 2006; Kashikaw a et al.l 
l2006bD . Deep broadband BVRi'z' and narrowband 
NB816 (Ac = 8I6OA, AAfwhm =120A) and NB921 
(Ac = 9196A, AAfwhm =132A) filter images were taken 
by the SDF project. All the images were convolved to 
have common seeing size of 0."98. Limiting magnitudes 
in 2" aperture at 3cr are (B, i?, i', z', NB816, NB921) = 
(28.45, 27.74, 27.80, 27.43, 26.62, 26.63, 26.54). Transmis- 
sions of these filters including CCD quantum efficiency, 
reflection ratio of the telescope prime mirror, the correc- 
tion for the prime focus optics and transmission to the 
atmosphere (airmass secz — 1.2) are also shown in Figure 

m 

Our NB973 image of the SDF was taken with Suprime- 
Cam mounted on the Subaru Telescope on 16 and 17 
March 2005. These two nights were photometric with 
good seeing of ~ 0."5-0."8. The total integration time 
is 15 hours. We have redu ced NB973 i mage fra mes us- 
ing the software SDFRED ^Ouchi et all [2004, : .Yagi et al.l 

Available from |http: / /soaps. naoj .org/ sdf / data/ 1 

IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 



I2OO2D in the same manner as m iKashikawa et all (|2004[ ). 
The NB973 image frames were dithered in a similar way as 
the SDF project did for other wavebands when they were 
taken. The combined NB973 image removed the slight 
fringing caused by OH-airglow that appeared in some im- 
age frames. The loss of survey area due to this dithering 
is only ~5%. The seeing size of the combined image was 
0."78 and convolved to 0."98, which is the common seeing 
size of the images of other wavebands, for the purpose of 

photometry^ S pectro photometric standard stars Feige34 

and Hz44 ((Okelll990D were imaged during the observation 
to calibrate the photometric zeropint of the stacked image, 
which is NB973= 32.03. The limiting magnitude reached 
NB973< 24.9 at 5cr with 15 hour integration. 

3. PHOTOMETRIC ANALYSIS 

3.1. Photometry 

After obtaining the stacked NB973 image, we conducted 
photometric analysis, making an object catalog. Source 
detection and photometry were carried out with S Extrac- 
tor software version 2.2.2 (jBertin fc Arnoutslll996l ). Pixel 
size of the Suprime-Cam CCDs is 0."202 pixel'^ We 
considered an area larger than contiguous 5 pixels with 
a flux [mag arcsec"^] greater than 2a to be an object. 
Object detection was first made in the NB973 image and 
then photometry was done in the images of other wave- 
bands using the double-imaging mode. 2" diameter aper- 
ture magnitudes of detected objects were measured with 
MAG_APER parameter while total magunitudes with 
MAC- AUTO. Low quality regions of CCDs, bright stellar 
halos, saturated CCD blooming, and pixels of spiky ab- 
normally high or low flux counts were masked in the SDF 
images of all wavebands, usi ng the official program 

cod«E3 

provided by the SDF team (jKashikawa et al.ll2004) . The 
final effective area of the SDF image is 876 arcmin^. The 
comoving distance along the line-of-sight corresponding to 
the redshift range 6.94 < z < 7.11 for LAEs covered by 
NB973 filter was 58 Mpc. Therefore, we have surveyed a 
total of 3.2 X 10^ Mpc^ volume using NB973 image. Then, 
the final object catalog was constructed, detecting 41,533 
objects down to NB973< 24.9 {5a). 

3.2. The Detection Completeness 

To understand how reliable our source detections are 
down to the limiting magnitude of NB973 < 24.9, we mea- 
sured the detection completeness of our photometry with 
the NB973 image. First, all the objects that satisfy our 
source detection criterion were removed from the NB973 
image using the SExtractor. Then, the starlist task in 
the artdata package of IRAlQ was used to create a sam- 
ple starlist of about 20,000 artificial objects with a random 
but uniform spatial and luminosity distributions ranging 
from NB973 — 20 to 25 mag. Next, using the mkobject 
task of IRAF, these artificial objects were spread over the 
NB973 image, avoiding the masked regions of the SDF and 
the locations close to the previously removed real objects 
with the distance shorter than 3/2 of their FWHMs. Af- 
ter this, SExtractor was run for the source detection in 
exactly the same way as our actual photometry. Finally, 
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we calculated the ratios of the number of detected artifi- 
cial objects to that of created ones to obtain the detection 
completeness. We repeated this procedure five times and 
averaged the obtained completeness. The result is shown 
in Figure [51 The completeness at our detection limit of 
NB973 = 24.9 is ^ 76%. The completeness was corrected 
when the number and luminosity densities of z = 7 LAEs 
were calculated in the §??. 

We evaluated the c ompl eteness in the same wa y as 
IShimasaku et all (|2006f) and iKashikawa et all (|2006bD did 
for the detection completeness of z = 5.7 and z = 6.6 
LAEs to keep consistency. However, in real life, some 
z = 7 galaxies will lie behind brighter sources at lower 
redshifts and thus the completeness correction will be arti- 
ficially small. The fractions that do can be included in the 
completeness estimate by simply adding artificial sources 
to the original image without masking out anything and 
without any exclusion zones in the placement of artificial 
sources. We also calculated our detection completeness in 
this way to see how much it is different from the original 
completeness evaluation. As expected, the completeness 
we calculated this time is slightly smaller. However, the 
difference is a factor of only 1.1-1.3 over NB973 — 20- 
25 and it does not change the evaluation of LAE number 
and luminosity densities much. Hence, for our consequent 
analyses, we use the original completesess calculated in the 
same way as done for z = 5.7 and z = 6.5 LAEs to keep 
consistency. 

3.3. Colors and Selection Criteria of z — 7 LAEs 

To isolate z — 7 LAEs from other objects, we investi- 
gated their expected colors and derived candidate selection 
criteria. We generated model spectra of LAEs at the red- 
shift ranging from z = 5 to 7 with rest frame Lya line 
equivalent width EWo(Lya) varying from to 300A as 
follows. First, we created a spectral energy distribution 
(SED) of a starburst galax y using a stellar populatio n syn- 
thesis model, GALAXEV (|Bruzual fc Chariot] [200l with 
a metallicity oi Z — Zq = 0.02, an age of t = 1 Gyr, 
Salpeter initial mass function with lower and upper mass 
cutoffs of tol =0.1 Mq and mu = 100 Mq and exponen- 
tially decaying star formation history for r = 1 Gyr. These 
parameters were chosen to be the same as those used to 
generate model z = 6.6 LAEs in iTaniguchi et al.l (|2005[ ) 
to keep consistency. Although recent observational stud- 
ies show that LAEs seem to be much younger than the 
1 Gyr age and/or 1 Gyr star formation d ecay time and 
z = 4.5 LAEs seem to have dust extinction (Gawiser et al.l 



I2006t iPirzkal et al.ll2007l : iFinkelstein et al.l 



2003), we did 



not consider the effects of dust on the SED since the two 
issues have opposite effects on the broadband colors of the 
LAEs and this does not have a major effect on the LAE 
selection criteria. Then, the SED was redshifted to each 
oi z = 5.0, 5.5, 5.7, and 7.0, and Lya absorp tion by 
IGM was applied to it, using the prescription of iMadaul 
(|1995() . Finally, flux of a Lya emission line with either of 
EWo(Lya) = 0, 10, 20, 50, 100, 150, 200, 250 or 300A was 
added to the SED at (1 + 2)1216A. We did not assume any 
specific line profile or velocity dispersion of Lya emission. 
Instead, we simply added 1/2 of the total line flux value, 
assuming that the blue half of the Lya line is absorbed by 
IGM. An example of a model spectrum of a 2 = 7 LAE is 



shown in Figure [T] 

Colors of these model LAEs were calculated using their 
SEDs and transmission curves of Suprime-Cam broadband 
and NB973 filters and plotted in a two-color diagram of 
z'— NB973 vs. i' — z' shown in Figure [31 As clearly 
seen in the diagram, a z = 7 LAE is expected to pro- 
duce significant flux excess in NB973 against z'. How- 
ever, it should be also noted that NB973 bandpass over- 
laps with the wavelength range at the longward edge of 
z' band. This allows LAEs and LBGs at even lower red- 
shifts z = 6.2-6.8 to cause the NB973 flux excess with 
respect to z' band if such galaxies have bright and steep 
UV continua. Actually, such objects were detected in 
our photometry. Their images and photometric proper- 
ties are shown and described in Figure |3| and t j3.4l Out 
of these lower redshift galaxies, z = 6.5-6.6 LAEs can 
be removed by requiring no detection in the narrowband 
fllter NB921 image whose bandpass corresponds to the 
Lya emissions at this redshift range iKodair a et 311120031 : 
ITaniguchi et al.|[2005l : IKashikawa et a l. 2006b). Hence, we 
classified NB973-excess objects with NB973 < 24.9 (5cr, 
2" aperture), which include our target z = 7 LAE, into 
following two categories based on z'— NB973 color. 

(1) z = 6.9-7.1 LAEs : B,V,R,i', NB816, z', NB921 
< 3a 

(2) z = 6.7-7.1 LBGs : B,V,R, NB816, NB921 < 3a, 
i' - z' > 1.3, z'- NB973 > 1.0 

where B, V, R, i' , NB816, z' and NB921 fluxes were mea- 
sured in total magnitudes while i' — z' and z'— NB973 col- 
ors in 2" aperture magnitudes. All the i' and z' aperture 
magnitudes fainter than 27.87 and 27.06 (2cr limits), re- 
spectively, were replaced by these values in the applica- 
tion of criterion (2). Since the flux of a LAE shortward 
of Lya emission should be absorbed by IGM, no detec- 
tions (< 3a) in B,V,R, NB816 and NB921 with either 
red i' — z' > 1.3 color or no detections in i' and z' were 
imposed as a part of the criteria. This can help eliminate 
interlopers such as L/M/T type dwarf stars and lower red- 
shift galaxies with other type of emission lines (e.g., H/3, 
[OIII], [OH], Ha, [SII] and so on). Also, criterion (1) im- 
plies that the robust z = 7 LAE candidates should show 
significant excess in NB973 over z' and NB921, z' ~ NB973 
> 1.72 and NB921 - NB973 > 1.64. 

Note that the color selection criteria (1) and (2) are 
slightly different from those in live et all (|2006( ) in that 
this time we include null detections in NB816 and NB921 
whose bandpasses correspond to Lya emission at z = 5.65- 
5.75 and 6.5-6.6, respectively, to make the criteria more 
reliable and sec ure. In fact, the object IOK-3 detected by 
live et all (j2006l) satisfied the criteria (1) and (2) simulta- 
neously except for NB921 < 3a and was spectrq s copical ly 
identified as a z = 6.6 LAE bv IKashikawa et al.l (|2006b[ ). 

This time, we found only one object s atisfying criterio n 
(1) (hereafter referred to as IOK-1 as in live et all (|2006[ )) 
and none met criterion (2). In order not to miss faint 
and diffuse z = 7 LAEs such as Lya blobs having ex- 
tended shapes with fairly bright cores but NB973 > 24.9 
(2" aperture mag.), we also loosened our detection limit 
cutoff adopting NB973 < 24.9 (total mag.) as another 
limiting magnitude. This increased the number of objects 
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satisfying criterion (1) by 17 while still no objects fell into 

However, this sample might be contaminated by some 
spurious objects such as sky residuals and noise due to 
fringing that might not be removed perfectly at the time 
of image reduction. Hence, we visually inspected all the 
broadband and narrowband images of each color-selected 
object and only kept those appearing to have condensed 
and relatively bright cores and excluded those having only 
diffuse faint shapes with no cores. More specifically, we 
removed objects that look apparently artificial such as 
connected bad pixels, tails of saturated pixels from bright 
stars and noises of discrete dismembered shapes or pieces 
of disconnected pixels with fairly large fluxes. As a result, 
we were l e ft wit h one object (hereafter called IOK-2 as in 
live et all ()2006[) ). The images of IOK-1 and -2 and their 
photometric properties are shown in Figure [¥] and Table 
[1] respectively. The color-magnitude diagram (z'— NB973 
vs. NB973) of IOK-1 to -2 as weU as afl the objects de- 
tected down to NB973 = 24.5 (total mag.) is plotted in 
Figure O Their two-color diagram is also shown in Figure 
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3.4. Possibility of Objects with Weak NB973-excess 
Being z ~ 7 LAEs 

As lTaniguchi et all (|2005f) did in selecting out their can- 
didate z = 6.6 LAEs in order not to miss faint targets, we 
also investigated the possibility of objects with a weak ex- 
cess of 1.0 > z'- NB973 > 3cr being z = 7 LAEs even 
though such objects do not have the expected colors of 
z = 7 LAEs predicted by the stellar populationin syn- 
thesis model in i i3.3l and Figure O We define the color 
criterion of such weak NB973-excess objects as: 

(3) B,V,R< 3a, i' - z' > 1.3, 1.0 > z'- NB973 > 3ct 

As mentioned in i )3.3[ our NB973 is located at the red 
edge of z' band. This could cause the criterion (3) to pick 
up interlopers such as z = 6.2-6.8 LAEs/LBGs (as Fig- 
ure [3] predicts), z = 1-3 extremely red objects (EROs) 
whose continua have the rest frame 4000A Balmer breaks 
that result in the NB973-excess against z' or M/L/T type 
red cool dwarf stars whose SEDs can have steep slopes at 
around NB973 bandpass. Such objects should be distin- 
guished from z = 7 LAEs if they are detected in NB973 as 
well as z'. From the photometry alone, it is difficult to tell 
if the criterion (3) objects are EROs or dwarfs. However, 
it is possible to say whether the objects are galaxies at 
z = 7 or not, which is more important in our study. 

According to the predicted colors of model galaxies in 
Figure [31 the criterion (3) should select out z — 6.2-6.8 
LAEs/LBGs, not z = 7 ones. However , as some of z = 
5.7 LA Es spectroscopically identified bv IShimasaku et al.l 
()2006[ ). though not so many, do not satisfy their color se- 
lection criteria computed using SED models, objects sat- 
isfying our criterion (3), which reside near the border of 
criteria (1) and (2), could be z = 7 LAEs. We found two 
objects to fall into criterion (3) (hereafter, referred as to 
Obj-4, the brighter of the two in 2" aperture NB973 mag., 
and Obj-5). Their colors, images and photometric prop- 
erties are shown in Figures [31 IH O and Table [H 

If they are LAEs, their redshifts can be further con- 
strained by using NB816 and NB921 images. As seen in 



Figure H Obj-4 is detected in i' , NB816, z', and NB921 
as well as NB973 but does not show any significant excess 
in NB816 against i' and in NB921 with respect to z' al- 
though it displays NB973-excess greater than 3(T against 
z'. Therefore, it is neither a z = 5.65-5.75 LAE nor 
z — 6.5-6.6 one. Since it is clearly detected in NB816, 
which is the waveband well shortward of z = 6.7-7 Lya 
emission, Obj-4 could be a z = 6.2-6.4 LAE or LBG. 

On the other hand, Obj-5 is detected in z' and NB921 
as well as NB973 but does not show significant excess in 
NB921 with respect to z' and thus is not a z = 6.5-6.6 
LAE. Also, detection in NB921 rules out the possibility 
of z = 6.7-7 LAEs since their fluxes shortward of Lya 
should be close to zero. Though it displays an excess of 
z'— NB973 > 3(7, its i' — z' color is very similar to that 
of a z ~ 5.7 LAE, which is predicted not to produce any 
NB973-excess. However, it is not detected in NB816 im- 
age and thus not a z ~ 5.7 LAE. Hence, Obj-5 could be a 
LAE or LBG at z = 6.2-6.4. 

As mentioned earlier, Obj-4 and -5 can be EROs or 
dwarfs. However, we have confirmed that all the objects 
with weak (> 3a) NB973-excesses (i.e., Obj-4 and -5) can- 
not be z = 7 LAEs and thus do not have to care anymore 
about the possibility of missing any faint z = 7 LAE can- 
didates. 

3.5. Possibility of IOK-1 and IOK-2 Being Variable 

Objects 

As the selection criterion (1) derived in H3.3\ shows, the 
most probable z = 7 LAE candidates are imaged in only 
NB973 waveband and not detected in any of other fil- 
ters. Since the NB973 image was taken 1-2 years after 
the BVRi'z' images of the SDF had been obtained, the 
sources only bright in NB973 can be some variable objects 
such as supernovea and active galactic neuclei (AGNs) 
that accidentally increased their luminosities during our 
NB973 imaging observation. Therefore, we investigated 
how many objects can be such variables. In another word, 
this corresponds to the number of the objects that were 
fainter than our detection limit NB973 = 24.9 (5cr) at 
some epoch but that can become brigher than it at an- 
other epoch. Since there are no enough data in z' and 
NB973 bands for the statistic of variables, we instead 
used i' band images ta ken over several separate epochs 
(jMorokuma et al.l l2007D for the best possible (but some- 
what rough) estimation we can do. 

First, we calculated the mean color of i'— NB973 over 
the range of NB973 = 22-25, which is < i'- NB973>= 
0.33, for the purpose of rough conversion of NB973 into 
i' magnitude. Using it, NB973 = 24.9 corresponds to 
i' = 0.33 -h 24.9 = 25.23. Since the detection limit of 
SDF i' band image (i' = 26.85 in 5a) is firmly deeper, the 
number count of objects fainter than our NB973 detection 
limit corresponding to i' — 25.23 can be securely obtained 
down to i' — 26.85. The number count per 0.5 mag bin as 
well as the magnitude increments needed to exceed NB973 
= 24.9 in brightness to be detected in NB973 are shown 
in Table [2] Since we were extrapolating the object num- 
ber counts in NB973 down to NB973 = 26.5 using i' band 
object number counts down to i' — 26.85 (5cr) and < i' — 
NB973>= 0.33, we also checked how similar the number 
counts in NB973 and in i' — 0.33 are to each other as shown 
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in Figure[S]and Table[3] Since i' — 0.33 number counts are 
slightly larger (by a factor of x 1.1-1.2 per bin), our cal- 
culation of the number of variables can be only a little 
overestimation. Note that the detection completeness of i' 
and NB973 are not corrected in their number counts. This 
can be the cause of the smaller counts in NB973 than i' 
toward our detection limit NB973 = 24.9. 

We use in our calculation four i' images of a part of 
SDF (~ 71% of the total area) taken at four separate 
epochs: 4 March 2005, 30 April 2003, 11 April 200 2 and 
24 April 2001, respectively (|Morokuma et al.ll2007D . The 
numbers of variable objects Ny{Ai') that increased their i' 
magnitudes by Ai' over the periods 2003-2005 and 2001- 
2002 were counted in each magnitude Ai' bin (matched 
to Am bin in Table ^ as shown in Table H In the 2003 
and 2005 images, ~ 70, 000 and 80,000 objects were de- 
tected down to their limiting magnitudes i' — 26.3 and 
26.6 (5cr, 2" aperture), respectively. Similarly, in the 2001 
and 2002 images, ~ 50,000 and 70,000 objects were de- 
tected down to i' = 25.9 and 26.2 (also, 5a, 2" aperture), 
respectively. Thus, taking the averages, we roughly as- 
sumed that Nobs = 75, 000 and 60,000 objects were de- 
tected in 2003-2005 and 2001-2002, respectively and di- 
vided the number of variables Ny{Ai') by these numbers 
Nobs to obtain the probabilities P(Ai') of finding the vari- 
ables with a brightness increase of Am — Ai' in the SDF 
down to our detection limit. 

Finally, multiplying the probability by the number 
counts of I'-detected objects N{Am) in Tableland sum- 
ming all them up, the number of variables that became 
brighter than NB973 = 24.9 came out to be 9-10. Note 
that since the magnitude increse of < Am < 0.1 is really 
small change and cannot be distinguished from photomet- 
ric errors in NB973 and i' magnitude measurements, which 
is also in the order of up to ^ 0. 1 , we ignored the number of 
variables in Am = 0-0.1 bin at the time of the summation. 
So far, we have considered only the data of the variables 
that increased their magnitudes over the two epochs and 
did not treat those having decreased their magnitudes. If 
we roughly assume that their numbers are approximately 
the same, the number of possible variables could be about 
one half of that we estimated above, which is ~ 4.5-5. 
Again, this number might be a little overestimation by a 
factor of 1.1-1.2 since we for our extrapolation used i' band 
number count instead of NB973 one, which is smaller as 
seen in Figure [5] and Table [31 Correcting for this factor, 
we estimate that the number of variables would be ^ 3.8- 
4. 5. This number is slightly different from that reported 
m live et all (|2006D since more elaborate calculations were 
used here. The estimated number of variables indicates 
that we cannot completely reject the possibility of nar- 
rowband excess of IOK-1 and IOK-2 being due to object 
variability. To securely reveal their identities, follow-up 
spectroscopy of them is required. 

4. SPECTROSCOPY 

To confirm the reality of our candidate LAEs, IOK-1 
and IOK-2, selected by the color selection criteria in i|3.3[ 
we carried out optical spectroscopy of them during 2005- 
2006 using the Faint O bject Camera And Spectrograph 
(|Kashikawa et al.l |2002l, FOCAS) on Subaru. The obser- 
vation status is summarized in Table [51 An Echelle grism 



(175 lines mm"^, resolution ~ 1600) with z' filter and 
0."8 slit was used to obtain the spectra of 30-min expo- 
sure each, dithered along the slit by ±1". The spectrum 
of spectrophotometric standard, either of Feige 34 Feige 
110 or BD-H28°4211 (lOkelll990l : [Hamuv et al.llT99l . was 
also obtained for each night and used for flux calibration. 
The observation data reduction a nd analysi s were all per- 
formed in the same manners as in live et all (I2006D . 

4.1. IOK-1, a z = 6.96 Lya emitter 

We identified IOK-1, the brighter of the two z = 1 LAE 
candidates, as a z = 6.96 LAE. The details of the spec- 
trosco pic analysis of this object were reported in live et al] 
(|2006l ). We measured the skewness and weighted skewness 
of the Lya emission line in IOK-1 spectrum and obtained 
S = 0.558 lb 0.023 and = 9. 46 ± 0.39 A, respectively . 
See lShimasaku et all (|2006D and lKashikawa et al.l (|2006bD 
for the definition of S and Sw These values show that 
the line is quite asymmetric and ensure that it is a Lya 
emission. Actually, our Sw value for IOK-1 is comparable 
to the average weighted skew ness of z = 5.7 and 6.6 L AEs 
(calculated from the data in IShimasaku et al.l (|2006t ) and 
iKashikawa erall (|2006bD ). < S^=^-'^ >= 7 A3 ± 1.47A and 
< 5"^=^ *^ >= 7.31 ± 1.5lA, respectively. 

The Lya line flux, i^(Lya), Lya line luminosity, L(Lya), 
the corresponding star formation rate, SFR{Lya) as well 
as other spectroscopic properties of the Lya emission 
line of IOK-1 are summarized in Table [6l To estimate 
SFR{Lya), we use t he following rela tion derived from 
Kennicutt's equation (iKennicutt 19981) w ith the case B 
recombination theory (|Brocklehurstlll97l[ ). 

S'i^i?(Lya) = 9.1 x lQ-'^^L{Lya)MQyr-^ (1) 
In addition, we estimate the UV continuum flux, 
F(\JY), by simply subtracting Lya emission line flux, 
i^(Lya), measured in the spectrum from NB973 to- 
tal flux, -FNBgys, obtained by SExtractor photometry 
(MAG- AUTO). 

F(UV) = i^NB973 - i^(Lya) (2) 
Then, this UV continuum flux can be converted into 
the UV continuum luminosity, Lj^(UV), and corre- 
sponding star formation rate, SFRjlJ V). To e stimate 
S F R(\]Y) , we use the following relation (IKennicu tt 1998; 
iMadau et aLlll998[ ). 

SFR{\JY) = 1.4 X 10~2^L^(UV)MQyr-i (3) 
The spectroscopic properties of the UV continuum of lOK- 
1 are listed in Table [3 

4.2. IOK-2 



As reported in live et all (|2006[ ). although there appears 
to be an extremely weak emission-like flux at around 
9750A (z = 7.02 if this is a Lya line) within the small gap 
between OH sky lines, 3 hours integration on the lOK- 
2 spectroscopy (obtained from 4 May 2005 and 24 April 
2006) was not deep enough to confirm if it is real or spu- 
rious since we had only S/N 2 even though measured 
within the gap. This did not allow us to draw any firm 
conclusion about IOK-2. 

To reveal the true entity of this object, we made 
additional 8 hours follow-up spectroscopy with Sub- 
aru/FOCAS on 10 April 2007 (See Table [S]). The seeing 
during this observing run was 0."4-l" with clear sky. We 
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combined the spectra taken at this night with those ob- 
tained in 2006 and 2005 to achieve the total of 11 hours 
integration. However, sky-subtracted stacked spectrum 
has shown neither the emission- hke flux at 9750A nor any 
other spectral features. We also combined only the spec- 
tra taken in 2007 and again could not find any emission 
lines. There are no signals that follow the dithering shifts 
among 30-min spectrum frames. This result indicates that 
the extremely weak emission- like flux at 9750A seen m 6 
hours stacked spectrum made from 2005 and 2006 frames 
is spurious. 

To see if our 11 hours spectroscopy has reached the 
depth required to detect a Lya emission, we compare the 
sky background RMS of the stacked spectrum with the 
Lya flux calculated from NB973 magnitude of the IOK-2. 
If we assume all of the flux in NB973 comes from the Lya 
line at z = 7 and adopt the total magnitude of NB973 
= 24.74 rather than 2" aperture one, we obtain the line 
flux of FP^°\Lya) = 2.9 x lO'^^ erg s'^ cm'^. On the 
other hand, binning of 4 pixels (corresponding to 0.017Mpc 
at z = 7) in the spatial direction is used to extract the 
one dimensional spectrum. The sky RMS (in terms of 
flux density) is measured in this spectrum by calculat- 
ing the variance in unbinned pixels along the dispersion 
direction within the wavelength range corresponding to 
NB973 passband 9655-9855A, and it is 3.0 x 10"^^ erg 
s~^ cm^^ A^^. The FWHM of Lya lin e, for example, of 
z = 6 . 6 LAE varies from 5.5 to 14. 6A (jKashikawa et al.] 
I2006bl : iTaniguchi et"ani2005l ). If we assume the FWHM 
distribution of z = 7 LAE is similar, then we obtain the 
Lya line flux of i^''P<''^(Lya) = (1.7-4.4) x lO'^^ erg s"! 
cm"^. This is 6.6-17 x fainter than FP^°*(Lya), indicat- 
ing that we have reached enough depth to detect the Lya 
line if IOK-2 is a real LAE at z = 7. Likewise, even if we 
use 2" aperture magnitude of NB973 = 25.51, we obtain 
FP'^°*(Lya) = 1.4 x lO'^^ erg s'^ cm'^, and F^'P^^Lya) 
is 3.2-8.2 X fainter than this. Furthermore, even if we as- 
sume ^ 68% of the NB973 flux comes from Lya line as 
live et all (|2006l) did, i^''P'='=(Lya) is still 4.5-12 (2.2-5.6) x 
fainter than i^P^°*(Lya) if we use NB973 total (2" aper- 
ture) magnitude to calculate i^P^°*(Lya). In all cases we 
have considered, our spectroscopy reached enough depth 
to detect a Lya line. Hence, IOK-2 might not be a LAE. 
However, we should note that residuals of the subtracted 
OH skylines around 9790 A in the 11 hours stacked spec- 
trum is still locally strong (~ 13% of the NB973 pass band 
is contaminated) and a Lya line can be masked out if it is 
weak and redshifted there. 

If IOK-2 is not a LAE, the possible origin of the NB973 
flux excess can be one of either a LBG at z ^ 7, a low-z 
ERO, a late- type star, a variable object or a noise. In the 
former three cases, spectroscopy could show no signals in 
the spectrum if their continuum light is very faint. If lOK- 
2 is a variable object, well possible as discussed in section 
13. 5[ it could have been fainter than our detection limit at 
the time of the follow-up spectroscopy. The possibility of 
IOK-2 being a noise spike in NB973 imag e cannot be rule d 
out though it is very low as described in live et al] ()2006f ). 
An additional NB973 imaging of SDF will be helpful to 
see if IOK-2 is either of a variable object or a noise. For 
the statistical study in the following sections, we hereafter 
consider that only IOK-1 is a z = 7 LAE we have success- 



fully identifled and IOK-2 is not. 

5. IMPLICATIONS FOR THE REIONIZATION AND GALAXY 
EVOLUTION 

From z ~ 6 quasar GP daigonostics, the neutral IGM 
fraction at this redshift was estimated to be x^^-^ ^ 0.01- 
0.04 and thus the reionization is believed to have already 
completed at around this epoch (jFan et al.l [2006') . This 
result is also supported by the spectral modeling analy- 
sis of the currently mo s t dist ant GRB at z ~ 6.3 con- 
ducted by iTotani et al.l (|2006D , placing the constraint of 
< a;gJ'6-3 < 0.17-0.6. 

On the other hand, the observed Lya LFs of LAEs at 
z ~ 6 and higher redshifts can be used to probe the epoch 
of reionization. The Lya LF is expected to decline beyond 
z 6 due to a rapid change of neutral IGM and ioniza- 
tion states before and after the completion of the reioniza- 
tion (iHaiman fc Spaans1ll999l : iRhoads fc Malhotra 11200 ll : 



iHu et al! 2002[ ). While the Lya LFs have been observed 
not to ev olve at z = 3-5.7 | Ouchi et al."2003'; 'Aii ki et al] 
2003; Tra n et al.|[2004t Ivan B rcukclcn et al. 2005), it was 
recently found to dechne as L*^g g ~ ^*=5 7 x (0.4-0.6) 
from z = 5.7 to 6 . 6 in SD F suggesting that x^^-^ < 0.45 
(jKashikawa et al.l l2006bD . Furthermore, we also found 
that the number density of z = 7 LAEs is only 18-36% 
of the density at z = 6.6 (jive et al.l |2006[ ). This series 
of decrements in densities might reflect the completion of 
reionization at around z ~ 6, beyond which the fraction 
of the neutral IGM hydrogen could possibly increase and 
attenuate the Lya photons from LAEs. 

However, this interpretation was based on the assump- 
tion that there had been no evoluion of LAE population 
from z = 5.7 to 7. The recent photometric study of z ~ 6 
i-dropouts and z ~ 7-8 z-dropouts in the Hubble Ultra 
Deep Field (UDF) demonstrated that galaxy number den- 
sity decreases by a factor of ~ 0.1-0.2, suggesting the 
rapid evolution of luminous ga laxies between these epochs 
(jBouwens fc Illingworthll2006D . 

In the following discussion, we re-evaluate the compar- 
ison of our LAE number and Lya luminosity densities at 
7 with those at z = 5.7 and 6.6, using the most 



up-to-date SDF data fr om IShimasaku et al.l (|2006t ) and 
iKashikawa et alj (|2006br ). We also investigate the possi- 
bility of LAE galaxy evolution between z = 5.7 and 7 and 
the degree to which it contributed to the number density 
deficit between these epochs. 

5.1. The Evolution of Lya LF at z > 6 

Figure |7| compares Lya line LFs at z = 5.7, 6.6 and 7 de- 
rived fr om the latest SDF data ( that is. Shimasaku et ~al\ 
(200^, IKashikawa et al.l (|2006bf ). and live et all (|2006D V 
In addition. Figure |5| shows the LAE number densities, 
nLya, Lya line luminosity densities, PLya, and correspond- 
ing star formation rate densities, SFRDi^ya, at 2.3 < z < 
7 down to our detection limit Liimit(Lya) — 1 .0 x 10"^^ 
erg s~ ^ (converted from NB973 < 24.9 (5ct) as live et al.l 
12006) did). Phya and SFRDi^ya at z = 7 are calculated 
using Lya line luminosity estimated from the spectrum of 
IOK-1 and equation [TJ The number and luminosity densi- 
ties at z < 7 are obtained by integrating the best-fit Lya 
Schechter LFs ((Schechter ,,1976.1 of LAEs down to our de- 
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tection limit Liimit(LyQ;) as follows. 



^Lya 



(f>{L)dL 



(4) 



(5) 



ilimit 



Phya — 



(j>{L)LdL 



( •J-'*'*-0.16' 



(6) 

42.60+^10, -1-5) for Lya LFs of LAE s at z = 5.7 and 6.6 in 
the SDF (taken from Table 3 in Kashi kawa et al. (2006b)), 
respectively. 0* = (22.0±12.0, 1.7±0.2) x lO-^Mpc'^ and 
L* = (5.4 ± 1.7, 10.9 ± 3.3) x 10^^ erg s"^ with a = -1.6 
are quoted fo r 2.3 < z < 4 5 (Ivan Breukelen et al.ll2005[ ) 
and z 4.5 (jPawson et al.l 12007 ) Lya LFs, respectively. 



We adopt {\og{(j3*[Mpc'^]),\og{L*[eTg 
3.44t°?°,43.04t°:i2,-l.5) and 
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We also use 
and L* = (5.8 
a = -1.5 for Lya LFs of LAE s at z = 3.1, 3.7 a nd 5.7 in 
- 1.0 deg2 of the SXDS field (|Ouchi et al.ll2007D . pLya is 
converted to SFRDi^ya by using equation [T] 

The uncertainties in the number and luminosity den- 
sities at z = 5.7-7 LAEs in the SDF in Figure [7] and [S] 
(and likewise Figure [H] and [TO]) include cosmic variance 
and the Poissonian errors associated with small-number 
statistic. To estimate the cosmic variance at z = 5.7-7, 
we adopt a bias parameter b = 3.4 ± 1.8 derived from the 
sample of 515 z = 5.7 LAEs detected in 1.0 deg^ o f 
the SXDS field (iSekieuchi et al.l[200l [O^chi et al.l[2005h . 
which is 5x wider than SDF. Then, applying the dark 
matter halo variances (z,crDM) = (5.7, 0.063), (6.6, 0.053) 
and (7.0, 0.044) obtained using analytic cold dark mat- 
ter model (jSheth fc Tormen 1 ll999HSomerville et al.l[200l 
and our SDF survey volumes to 6 = CTu/ctdm, we calculate 
the geometric mean of cosmic variance at z = 5.7-7, which 
is 8.4-27%. The maximum cosmic variance of = 27% 
is included in the errors in Figure FTHTOl Similarly, the cos- 
mic variance at z = 2.3-4.5 is also calculated and included 
in the Figure 8. The Poissonian errors for smal l-number 
statistic are estimated using Table 1 and 2 in IGehrels 1 
((1986) . When the densities and errors are calculated for 
z = 5.7, 6.6 and 7 LAEs in the SDF, the detection com- 
pletenesses in NB816, NB921 and NB973 images are also 
corrected (see Figure [5] for NB973 completeness). 

While it remains unchanged at 2.3 < z < 5.7, the LAE 
number density decreases by a factor of "-LyQ ^/"•Lya '' — 
5.7 to 6.6 and n£=I/n£=^-S ~ 0.17 from 
Similarly, the LAE Lya luminosity den- 
"^^-^ I Pl=l-'' - 0.21 and 



0.24 from z = 
z = 6.6 to 7. 

sity dechnes by factors of pl^^ /^Lya 
Phya/ Phya ^ ~ 0.15. If we assumc that the LAE popula- 
tion does not evolve from z = 7 to 5.7, this density deficit 
might reflect an increase in neutral IGM hydrogen with 
redshifts. 

However, the density decline might also possibly be as- 
cribed to the evolution of LAE population. If the number 
of LAEs having luminosities fainter than our SDF detec- 
tion limits drastically increases from z = 5.7 to 7, this can 
certainly affect our estimations of nj^ya and pLya- Hence, 
the Lya LF alone cannot resolve this degeneracy between 
the reionization and galaxy evolution effects. 



To cope with this matter, the rest frame UV contin- 
uum luminosity function (UVLF) of LAEs can be used 
to extrtact the galaxy evolution effect alone si nce it is 
not sup pressed by neutral hydrogen. Kashik awa et al.l 
(|2006bli have compared the UVLF (rest frame ~ 1255A at 
z = 6.6 and - 1350A at z = 5.7) of LAEs in SDF and 
other field also imaged by Suprime-Cam and found that it 
does not significantly change from z = 5.7 to 6.6. This 
suggests that the density deficit between z = 5.7 and 
6.6 are not mainly cau sed by galaxy evolution. Thus, 
I Kashikawa et al.l (]2006tf ) concluded that the reionization 
might have ended at around 5.7 < z < 6.6 and it supports 
the results of z ^ 6 quasars and GRB ((Fan et al.i 120061: 
iTotani et aLll2006f) . If this is also the case for z = 6.6-7 
LAEs, the further decline of LAE density implies increase 
in nuetral hydrogen that attenuates Lya photons and sup- 
ports J£ashi^awajvt_al^ (2006bV s result. 

5.2. Can the Lya LF evolution be explained only by 
galaxy evolution? 

We do not know whether the LAEs themselves evolve 
at 2 = 6.6-7. If the galaxy evolution occurs at z = 6.6- 
7, the further decline of LAE density at these epochs re- 
flects it in addition to reionization. Hence, in this section, 
we investigate the possibilities of the LAE evolution from 
z = 6.6 to 7 using three independent methods: (1) Com- 
parison of the UVLFs of z = 5.7 and 6.6 LAEs with that of 
z = 7 LAEs derived from our spectroscopic data of lOK- 
1, (2) Estimation from the UVLF evolution of LBGs and 
(3) Application of an LA E evolution model constructed 
bv lKobavashi et aT] (|2007D b ased on a hierarchical clust er- 
ing galaxy formation model (|Nagashima fc Yoshiill200^ to 
predict the expected change of Lya LF from z = 7 to 5.7 
due to galaxy evolution alone. 

5.2.1. Implications from UVLF of z — 7 LAEs 

First, we roughly estimate UVLF of z = 7 LAEs to see 
if there is any possible galaxy evolution from z = 6.6. 
We calculate absolute UV magnitude Afuv,i230 at the 
rest frame 1230A for IOK-1 from the UV continuum flux 
F(UV) obtained in H4.1l using equationj^jand F(Lya) mea- 
sured in the spectrum of IOK-1. That is, 



-A^UV,1230 — "^UV,1230 ^ DAI 



-2.5 log 



cAA 



2.51og(l 
F(UV) 



6-L»M-|-2.51og(l-Ki) 



where touv,1230 is UV apparent magnitude, A — 1230(1 + 
z)A and AA is the wavelength range in which the UV con- 
tinuum is covered by NB973 passband, AA — 9855A — 
(1 -|- z)1216A. Also, DM and c are a distance modulus 
and the speed of light. Figure |9| shows the UVLF of z = 7 
LAE derived here together with those of z = 5.7 and 6.6 
LAEs. We ignore a subtle difference in the rest frame UV 
wavelengths (rest frame ^ 1230A at z = 7.0, ^ 1255A at 
z = 6.6 and ^ 1350A at z = 5.7) assuming that the LAEs 
have flat UV continua. Also, the detection completeness 
of NB973 image is corrected using Figure ID The UVLF 
implies that there is no galaxy evolution from z = 7 to 6.6, 
and the density deficits of riLya and pLya between these 
epochs might be attributed mainly to the reionization. 
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5.2.2. Estimation from the UVLF Evolution of LBGs 

Even though the z = 7 UVLF derived from the SDF 
data suggests that LAEs do not evolve from z — 7 to 6.6, 
it suffers small statistics due to the relatively shallower de- 
tection limit in NB973 (equivalent to L(Lya) > 1.0 x 10^'' 
erg s~^). Therefore, we discuss the possibilities of the LAE 
galaxy evolution at z = 5.7-7 using inferences from other 
independent methods, by which we try to obtain some 
helpful insights. 

One possible way to estimate the LAE galaxy evolu- 
tion at z = 5.7-7 is the inference from the evolution of 
UVLF of high-z LBGs, assuming LAEs and LBGs share 
a similar evolutionary history. We use the UVLF data 
from the recent observational studies abou t z ^ 4-8 LBGs 
con ducted bvlYoshida ct al (2006), Bou wens et al.l (|2006[ ) 
and lBouwens fc lUingworthl pOOet ). Their surveys, when 
combined together, form the deepest and the widest imag- 
ing data with the samples of the largest numbers in all 
LBG surveys. Interestingly, lYoshida et al.l (12006) com- 
bined their data of z ~ 4 and 5 LBGs with those from 
lo wer-z LBG surveys an d z ~ 6 LBG (i-dropout) study 
bv lBouwens et all |2006f ) and found clear evolution of the 
UVLF from z ~ 6 to 0, in which only the characteris- 
tic magnitude, M^y, changes significantly and almost lin- 
early with redshift while the normalization factor, 0*, and 
the faint end slope, a, tend t o rem ain constant as seen in 
Figure 22 of lYoshida eral] (|2006l ). This trend of M*^, 
(j)* and a vs. z continues up to z ^ 7.4 when we add 



I's Lya line luminosity, yield 77,°'^?°'^*'^ 



M^y = -19.5 ±0.6 mag 



= 0.00202t° !^™?^ Mpc-3 and 



-1.73 (or M, 



uv 



00076 

-18.75 ± 0.6 mag. 



= 0.00218 



Mpc-3 and a = -1.7 3) of the first (or second) LBG UVLF 

ling wort hi ([1 ' 
between z 



at z ~ 7.4 derived by iBouwens &: Illingwort (pool . 
We estimate the change in M^y between z = 5.7 and 

4, 5, 6 



z = 7 from the z-dependence of M/jy at z 



and 7.4, assuming the correlation is linear with the slope 
of AMyy/Az ~ 0.47 mag. As a result, Myy is ex- 
pected to become fainter by 0.6 mag, which corresponds 
to the luminosity of L^f^"^* ^ L*^5.7 x I0-0-4xo.6 ^ 
-^z=5 7 X 0.58. Here, the relation between the equations 
H]and the Schechter LF in absolute magnitude form. 



<j){M)dM = -<?!)*(lnlO) 
5 



is used to convert A^uy to L 



exp 



.lOi(M'-M) 

(8) 

To infer the deficit by 
which Lya LF of LAEs decreases from z = 5.7 to 7 due 
to thier evolution alone, we now roughly assume that 
this UVLF evolution of the LBGs can be also applied 
to that of LAEs at z = 5.7-7 and Lya line luminosities 
of LAEs are simply proporti onal to their UV conti nuum 
luminosities as Figure 15 in iTaniguchi et"al] (|2005[ ) sug- 
gests. Based on this idea, we change i*=5.7 of our best-fit 
Schechter Lya LF at z = 5.7 in exactly the same way (i.e., 
log((^*[Mpc-3]) = -3.44t° jo and a = -1.5 as in ^but 
log(L*^5 7[erg s-i]) = 43.04^°;}^ + log 0.58 this time) to 
obtain z = 7 Lya LF. This result is compared in Figure 
[TO] with actual observation data of IOK-1. 

The inferred Lya LF at z = 7 does not really agree 
with one calculated from the spectrum of IOK-1. Our 
density deficit between z = 5.7 and 7 LAEs cannot be 
explained by only the galaxy evolution factor estimated 
here. The integrations of the inferred Lya LF using equa- 
tion [5] and [5] down to logL(Lya) = 43.05, which is lOK- 



expect,^— 7 



'Lya 

2.3 X 10^* 



erg s' 



i 1.5 X 10~5 
"iMpc-3, 



re- 



Mpc-^ and p^/J 
spectively. Our LAE number and Lya line luminosity 
densities at z = 7 are nf=^ ~ '^^ ^ m-e ^^,^-3 

and Pl^*" 



{A.ltll) X 1037 



'Lyo - (3-6+y) X 10-« Mpc- 
_ ^„ erg s~^Mpc~'^ based on lOK- 

1 data alone, respectively. Therefore, the density deficits 

r z=7 /„expect,z=7 ^ r)^+0A9 i z=7 / expect,z=7 ^ 
01 ^Lya/^-Lya — ^-^^-0.19 ^^'^ Phya/ Phya — 

0.181q;'^3 might be due to the attenuation of Lya photons 
by neutral IGM having existed during the reionization. 

In order for the inferred Lya LF at z = 7 to have the 
same number density as observed Lya LF at z = 7 (i.e., 



n 



expect, z— 7 

Lya 



Ly. 



'^), we have to change the characteristic 



luminosity 7^*^^?'='="^ a factor of x0.65lg'^g. This factor 
might reflect the deficit by which Lya LF of LAEs de- 
creases from z = 5.7 to 7 due to the attenuation of the 
Lya lines of LAEs by the increasing neutral IGM during 
the reionization beyond z ^ 6. We can refer to such a 
deficit factor due to the neutral IGM attenuation as IGM 
transmission to Lya photons, Tl^^. In the case of our dis- 
cussion so far, this can be regarded as the ratio of the Lya 
line luminosities of LAEs in the environments with some 
neutral IGM fraction xhi still remaining and with no neu- 
tral IGM (i.e., xm = 0), Tl™ = L'^^^ihya)/ L'=^^=^{hya). 
Once we know Tl™, the neutral IGM fract ion a± z = 7, 
be estimated. However, the calculation of xm 
from Tly^ is not a simple issue and is dependent on the- 
oretical models. We will discuss it in H5.3I 

Similarly, the z-dependence of M^y, AMyy/Az ~ 0.47 



mag, predicts AM, 



T expect 
^z=6.6 



UV 



,.7 X 10 



0.42 for Az = 6.6 - 5.7 and thus 
L*^5 7 X 0.68 due to 



-0.4x0.42 



LAE galaxy evo l ution f rom z = 6.6 to 5.7. However, 
iKashikawa et al.l (|2006b[ ) found that the Lya LF declines 
in such a way that Lz=6.6 ~ ^z=5.7 ^ (0.4-0.6) from z = 5.7 
to 6.6, regarding their photometric and spectroscopic LFs 
as the upper and lower limits of z = 6.6 Lya LF, respec- 
tively. Hence, the attenuation of Lya photons by the neu- 
tral IGM at z = 6.6 is Tl™ ^ 1%^^ JKtlT' - 0-59- 
0.88. 

^-j^j The decrease in Lya LF from z — 5.7 to 6.6 and 7 can- 
not be explained only by the evolution of LAEs inferred 
from that of LBGs. This result implies that the remain- 
ing deficits could come from the attenuation of Lya lines 
by neutral IGM. If this is the case, Lya line tends to be 
more attenuated at higher redshift as ~ 0.59-0.88 at 

z = 6.6 and 0.65^Q'^g at z = 7, implying that the neutral 
IGM fraction, xhi increases with redshift beyond z '--^ 6 
as derived in t J5.3l However, note that this result is based 
on the assumption that LAEs evolve in the same way as 
LBGs do. This might not be necessarily true. Although 
the LAEs arc believed to be closely related to LBGs and 
many of candidate LBGs at high redshift have been iden- 
tified as LAEs by spectroscopy, the link between these two 
populations has not been clearly understood yet and they 
might have followed different evolutionary histories. 

5.2.3. Application of A Galaxy Evolution Model 

In the previous sections, we tried to estimate the intrin- 
sic evolution of Lya LF of LAEs from UVLF evolution 
of LAEs and LBGs, with an implicit assumption that the 
evolutions of Lya and UV luminosities are similar to each 
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other. However, this assumption may not be true in real- 
ity, and hence our argument will be strengthened if we can 
show that these are indeed similar in a realistic theoretical 
model of LAEs. 

For this purpose, we use a recent mo d el for LAE evo- 
lution constructed by iKobavashi et al.l (j2007l . hereafter 
K07). This model is an extension of one of the lat- 
est hierarchical clustering models of galaxy formation 
(Na ga shim a fc Yoshii 2 004) , in which the merger histo- 
ries of dark matter haloes are modeled based on struc- 
ture formation theory and star formation processes in dark 
haloes are calculated to predict the photometric properties 
of galaxies. This model can reproduce most of the obser- 
vations for photometric, kinematic, structural, and chem- 
ical properties of local galaxies, as well as high-z LBGs 
(|Kashikawa et al]|2006a[ ). K07 extended this model with- 
out changing the original model parameters, but intro- 
ducing new modeling only for the escape fraction of Lya 
photons ifcZ") from galaxies based on physical considera- 
tions. Specifically, the dust extinction of Lya photons and 
effect of galaxy-scale outflow are newly taken into account. 
This is the first model for LAEs based on a hierarchical 
galaxy formation model in which /^^" is not a universal 
constant but depends on physical conditions of galaxies. 
This model can reproduce the observed Lya LF of LAEs 
in z ~ 3-6, and predicts that galaxies under strong galaxy 
scale outflow with f^'^." ^ 1 are dominant in the bright- 
end of Lya LFs, which is also consistent with observations. 
It should be noted here that /^^" in the K07 model could 
vary from galaxy to galaxy and may evolve within a galaxy, 
and hence even if the Lya photon production rate is pro- 
portional to star formation rate, the evolutions of Lya and 
UV LFs could be different from each other. 

The K07 predictions of the Lya and UV LFs of LAEs 
at z = 5.7, 6.6 and 7 assuming T^y^ = 1 are presented 
in Figure [7] and [SI respectively. The evolutions of number 
density and Lya luminosity density of LAEs with a thresh- 
old Lya luminosity predicted by this model are shown Fig- 
ure [H As demonstrated in K07, the deficit of the observed 
LAEs compared with the model prediction of Lya LF is 
clear at z > 6 as seen in Figure [7] while this model pre- 
cisely reproduces the observed evolution at z '^3-6. On 
the other hand, the degree of evolution of UVLF of LAEs 
predicted by the model is similar to that observed in the 
same redshift range. 

The fact that the model prediction is consistent with the 
UVLF evolution but not with the Lya LF evolution then 
implies that the evolution of the observed Lya LF at z < 6 
could be caused by the IGM absorption. The discrep- 
ancy can be resolved if we adopt a simple prescription of 
luminosity-independent IGM transmission: = 0.62- 

0.78 at z = 6.6 and Tl^^ = 0.40-0.64 at z = 7. 

5.3. Implications for Reionization 

In the previous section we have shown that the evolu- 
tion of Lya LF at z > 6 could be likely a result of Lya 
photon absorption by neutral IGM, implying a significant 
evolution of the IGM neutral fraction beyond z > 6. In 
order to obtain some quantitative implications for reion- 
ization, however, we must translate the estimates of T^^^ 
obtained in the previous section into the IGM neutral frac- 
tion, a;Hi- This procedure is not straightforward because 



this translation is generally model dependent (e.g., Santos 
2004; Dijkstra et al. 2007). 

Here, we apply the dynamic model with a reasonable 
velocity shift of Lya line by 3 60 km redward of the 
systemic velocity ()Santosll2004D . The attenuation factor of 
Lya luminosity is given as a function of xm, and the rea- 
son for the choice of this model is that this model predicts 
no attenuation when xm — 0. Note that some other mod- 
els of Santos (2004) predict a significant attenuation even 
in the case of xhi = 0, due to the neutral gas associated 
with the host haloes of LAEs. Choosing this particular 
model then means that we ascribe the evolution of the 
Lya LF at z > 6 only to the absorption by pure IGM. We 
consider that this is a reasonable assumption, since obser- 
vations indicate that the escape fraction o f Lya photons 
is about unity at least for LAEs at z 3 (Gawi ser et al.l 
2006). If LAEs at z 7 are a similar population to the 
low-z LAEs, we do not expect significant absorption by 
neutral gas physically associated to LAEs. On the other 
hand, it should also be kept in mind that if z ^ 7 LAEs 
are surrounded by a significant amount of nearby neutral 
gas that is not present for low-z LAEs, the estimate of xhi 
as an average of IGM in the universe could become lower 
than those derived here. 

In section [5X2l we obtained Tl'^^ = 0.59-0.88 and 

0- 65tn ?g at z = 6.6 and 7.0, respectively. Apphcation 
of iSantod ()2004[ ) model yields the neutral fractions of 
xl^^-^ ~ 0.12-0.42 and xgf/'j - 0.12-0.54. If we use 
the K07 model in section [5.2.31 to estimate xhi, we find 
xfiY^'^ - 0.24-0.36 from Tl'^^ = 0.62-0.78 at z = 6.6, and 
xfiY^ ~ 0.32-0.64 from Tl^^ = 0.40-0.64 at z = 7. The 
neutral fraction xm at z = 6.6 and 7 estimated from two 
independent methods are consistent with each other and 
tends to increase with redshift at z > 6. These series of 
a;Hi values at z = 6.6 and 7, combined with x^^ '^ ^ 0.01- 
0.04 and xf^^-^ < 0.17-0.6 derived from quasar GP tests 
and G RB spectral analysis f|Fan et al.l 120061 : iTotani et all 
[20061) . supports the picture that the reionization com- 
pleted at z ~ 6, beyond which it was still in progress with 
larger neutral fraction of IGM hydrogen, which evolved 
with redshift. The neutral IGM fractions obtained by in- 
dependent methods are summerized in Table [5] 

However, our constraint suggests that the neutral IGM 
persists at the ~ 50% level as late as z = 7 and would con- 
tradict the WMAP conclusion that the reionization epoch 
was z — 10.9^2 3 at the > 95% confidence level. Our re- 
sults could be reconciled with WMAP only if there is a 
statistical fluke (one time in 20, a 95% confidence range is 
wrong) or the reionization happened twice (e.g., Cen 2003) 
so that a lot of the observed electron scattering happens 
at z ^ 7, and then the universe becomes partially neutral 
again, allowing us to observe neutral gas at z = 7. 

Finally, we again emphasize that these quantitative re- 
sults are model-dependent and should be interpreted with 
caution. However, the decrease in the Lya LF of LAEs be- 
yond z ~ 6 is more significant than expected from UVLF 
evolution or a theoretical model, and hence the physical 
status of IGM might be changing at z > 6. 

6. SUMMARY AND CONCLUSION 

We have conducted a narrowband NB973 survey of z = 
7 LAEs, established color criteria to select out z = 7 LAEs, 
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and found two candidates down to L{Lya) > 1.0 x 10^^ erg 
s^^ {5a). By follow-up spectroscopy, the brighter of the 
two was indentified as a z = 6.96 LAE while we can con- 
firm neither Lya emission nor any other features in the 
spectrum of the other candidate despite the sufficiently 
long integration time. 

The number and Lya luminosity densities at z = 7 ob- 
tained by this study were compared to those at z = 5.7 and 
6.6 deriv ed from the latest samples obtained by the SDF 
surveys (jShimasaku et al.l |2006[ iKashikawa et all I2006H ) 
down to our detection limit, and clear evolution of den- 
sity deficits with increasing redshifts was observed such 
that: nl=;^J/nlf^-^ ~ 0.24 and n£;^/„£=6.6 ^ q^j. 

p17'-Vp17^-' ^ 0.21 and plfjpi^t-' ^ 0.15. If we as- 
sume that the LAE population does not evolve from z = 7 
to 5.7, this series of density deficits could reflect an increase 
in neutral IGM hydrogen with redshifts beyond z ~ 6. 

To see if LAE evolves from z — 7 to 6.6, we also com- 
pared UVLF of z = 7 LAE with those of z = 5.7 and 6.6 
LAEs derived from the SDF LAE surveys. No decrease in 
the number density of UVLF was observed from z = 5.7 
through 6.6 to 7. Since the UV photons are not attenuated 
by neutral IGM and the UVLF is only sensitive to galaxy 
evolution, our result suggests that the deficits in JiLya and 
Phya might reflect the cosmic reionization and the LAE 
population does not significantly evolve at z = 5.7-7. 

However, the UVLF at z = 7 suffers from small statistics 
at this time and the interpretation is not robust. Hence, 
the amount by which the LAE evolution affects the density 
deficits among z = 7, 6.6 and 5.7 wer e investigated by the 
inference from UVLFs of z < 8 LBGs (lYoshida et al.l l2006l: 
iBouwens fc IllingworthI l2006t iBouwens et al.ll2006D based 
on the assumption that LAEs would have evolved in the 
same way as LBGs and Lya line luminosities of LAEs are 
proportional to their UV continuum luminosities. Even af- 
ter the galaxy evolution was taken into account, there still 
remained some density deficits among these epochs. If we 
attribute the deficits to the attenuation of Lya photons by 
the neutral IGMs, the neutral fractions of the Universe at 
z = 6.6 and 7 are estimated to be 0.12-0.42 and 0.12-0.54, 
respectively. This result, combined with neutral fractions 
derived from z ~ 6 quasars and a z ~ 6.3 GRB, supports 
the completion of the reionization at z ~ 6 and the possi- 
ble evolution of neutral IGM beyond this redshift. 



Again, this result is based on the assumption that LAEs 
would have evolved in the same way as LBGs, which 
might not be necessarily true. Therefore, we furthermore 
used a LAE evolution model (K07 model) constructed 
from hierarchical clustering scenario to reproduce Lya LFs 
at z ~ 5.7, 6.6 and 7 in the case of transparent IGM 
[xm = 0) and com pared with Lya LFs obtained by the lat- 
est SD F survey s ( Shimasaku et al. 2006'; Kashika wa et al.l 
■2006bt live eTal., 200 6). The observed data at z = 6.6 and 
7 showed smaller number and luminosity densities than 
those predicted by the model, suggesting that there still 
remains the possibility of the incomplete reionization at 
those epochs. The neutral fractions at z = 6.6 and 7 es- 
timated from the decline of the LFs by the reionization 
factors alone after the galaxy evolution effects had been 
corrected are a;gf^-6 ~ 0.24-0.36 and xgf^ - 0.32-0.64, 
respectively, also consistent with quasar and GRB results. 

The results regarding z = 7 LAE presented here are 
based on relatively shallow depth of NB973 imaging, small 
sample statistics, and only the SDF and optical imaging 
data. From these data alone, the trend of the density 
deficit between z = 5.7 and 7 in fainter LAE populations 
and in other places in the Universe, changes in physical 
properties of LAEs associated with their evolution be- 
tween these epochs, and typical spectroscopic properties 
of z = 7 LAEs such as the direct detection of the attenua- 
tion of Lya line cannot be inferred. Also, the calculation 
of the UV continuum flux and thus UVLF is dependent on 
relatively rough estimation without infrared data. Deeper 
NB973 imaging of the SDF as well as other fields for which 
infrared images are available and follow-up spectroscopy of 
newly detected LAE candidates will provide the answers 
and more precise results in future studies. 
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We express the gratitude to the SDF team for obtaining 
and providing us with invaluable imaging data. K.O. ac- 
knowledges the fellowship support from the Japan Society 
for the Promotion of Science and the Special Postdoctoral 
Researchers Program at RIKEN. 
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Table 1 



Photometric properties of candidate Lyo emitters and NB973-excess objects 



Object 


Coordinate 


i' 


NB816 


z' 


NB921 


NB973 


NB973 (total) 


Criteria 


IOK-1'" 


13:23:59.8 +27:24:55.8 


>27.84 


>27.04 


>27.04 


>26.96 


24.60 


24.40 


(1) 


IOK-2 


13:25:32.9 +27:31:44.7 


>27.84 


>27.04 


>27.04 


>26.96 


25.51 


24.74 


(1) 


IOK-3'' 


13:24:10.8 +27:19:28.1 


>27.84 


>27.04 


26.26 


25.08 


24.87 


24.57 


C 


Obj-4'* 


13:25:09.1 +27:32:16.8 


27.14 


26.77 


25.75 


25.51 


24.97 


24.85 


(3) 


Obj-5'* 


13:23:45.8 +27:32:51.4 


>27.84 


>27.04 


25.76 


25.44 


25.10 


24.74 


(3) 



NOTE: Units of coordinate are hours: minutes: seconds (right ascension) and degrees: arcminutes: arcseconds (declination) using J2000.0 
equinox, i' , NB816, z' , NB921 and NB973 are all 2" aperture magnitudes while NB973 (total) is a total magnitude. Magnitudes are replaced 
by their 2a limits if they are fainter than the limits. Color criteria are those used to select out candidate LAEs and classify z'— NB973 excess 
objects (See ^331 and l3?it . 

"IOK-1 was proven to be a 2 = 6.96 LAE spectroscopically by live et al.l l|2006l ) . 

''IOK-3 was found to be a 2 = 6.6 LAE by independent spectroscopy bv lKashikawa et al.l Il2006bl) 

'^IOK-3 satisfies color criteria (1) and (2) simultaneously except for NB921 < 3a. 

"^These objects show flux excess of 1 > z'~ NB973 > 3a but could be either of 2 = 6.2-6.4 LAEs, low-2 ellipticals or latc-type stars (See i|3.4l l. 



14 



Ota et al. 



Table 2 

The number counts op ^'-detected objects fainter than our detection limit NB973 = 24.9 



i'" 


NB973" 


Am" 


N{i') = NiAmY 


25.23 


24.9 (5cr) 


0.0 




25.23-25.33 


24.9-25.0 


0.0-0.1 


2297 


25.33-25.83 


25.0-25.5 


0.1-0.6 


13589 


25.83-26.33 


25.5-26.0 


0.6-1.1 


16012 


26.33-26.83 


26.0-26.5 


1.1-1.6 


17367 



NOTE: Units of the first to third columns arc all AB magnitudes with 2" aperture, 
"j' band magnitudes were converted into NB973, using < i'— NB973 >= 0.33. 

''The NB973 magnitude increments required for objects to become brighter than our detection limit NB973 (Scr) = 24.9. That is, Am = NB973 
-24.9 or Am = i' - 25.23. 
°The number counts of i'-detected objects. 
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Table 3 

The number counts of objects in the NB973 and i! images of the SDF 



NB973 = 


%' - 0.33° 




iV(NB973)'= 


N(i'Y 


22.5 


23.0 


22.83-23.33 


2599 


2588 


23.0- 


-23.5 


23.33-23.83 


3808 


3953 


23.5- 


-24.0 


23.83-24.33 


5230 


5697 


24.0- 


-24.5 


24.33-24.83 


7093 


8121 


24.5- 


-25.0 


24.83-25.33 


8800 


10660 


25.0- 


-25.5 


25.33-25.83 




13589 


25.5- 


-26.0 


25.83-26.33 




16012 


26.0- 


-26.5 


26.33-26.83 




17367 



NOTE: Units of the first and second columns ajre total magnitudes, 
"i' band magnitudes in the second column were converted into NB973 magnitudes in the first column, using < i' — NB973 >= 0.33. 
^i' magnitudes brighter then -Scr limiting magnitude of the SDF i' band image i' (5cr) = 26.85. 
"^Number counts of objects detected in the SDF NB973 image down to our detection limit NB973 (5<t) = 24.9. 
'^Number counts of objects detected in the SDF i' image down to its limiting magnitude i'{5a) = 26.85. 
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Table 4 

The number counts of i'-DETECTED VARIABLES AGAINST i' BRIGHTNESS INCREMENT 





Number of 


variables Nr(Ai') 


P(Ai'Y' = 


Nr(An/N„,, 


P(Ai') > 


< N(Arn) 


(AB uirtg) 


2003 2005 


2001 2002 


2003 2005 


2001 2002 


2003 2003 


2001 2002 


0.0 














0.0-0.1 


250 


409 


3.3 X 10-'^ 


6.8 X 10--'' 


7.6^= 


15.6^= 


0.1-0.6 


52 


37 


6.9 X lO-'' 


6.2 X lO-'^ 


9.4 


8.4 


0.6-1.1 


1 


2 


1.3 X 10-5 


3.3 X 10-5 


0.21 


0.53 


1.1-1.6 


1 


1 


1.3 X IQ-s 


1.7 X 10-5 


0.23 


0.30 



Number of variables that became brighter than NB973 = 24.9 9.8'^ 9.2'^ 



NOTE; The number of variables brighter than our detection limit in the stacked i' SDF image, i' = 0.33+ NB973 {5a) = 25.23, were counted. 
"Increase in i' magnitude over the two epochs, which is binned to match the third column (Am) of Table 3. 
''Probability of finding a variable with a brightness increase of Ai' = Am in SDF down to our detection limit. 
°We ignore these values since the Am = 0-0.1 cannot be distinguished from photometric errors. 
'^These figures were obtained by EAm-P{Ai') x N{Am) = EajP (A?') x N{Am). 



z = l LAE 



Table 5 

Status of the follow-up spectroscopy 



Object 


date 
(HST) 


seeing 
(") 


exposure" 
(seconds) 


FOCAS Mask 


lOK-l" 


14, 15 May 2005 


0.5-0.7, 0.9-1.0 


10800 


MOS-1 




1 June 2005 


0.6-0.8 


3600 






24 April 2006 


0.9-1.5 


16200 




IOK-2^ 


14, 15 May 2005 


0.5-0.7, 0.9-1.0 


3600 


MOS-4 




24 April 2006 


0.9-1.1 


7200 






19, 21 June 2006 


1.0-2.0, 1.0-2.0 


28430"^ 






10 April 2007 


0.4-1.0 


28800 




IOK-3 


14, 15 May 2005 


0.5-0.7, 0.9-1.0 


3600 


MOS-2 




1 June 2005 


0.6-0.8 


5400 




Obj-4'= 


14, 15 May 2005 


0.5-0.7, 0.9 1.0 


1800 


MOS-5 


Obj-5^ 


14, 15 May 2005 


0.5-0.7, 0.9-1.0 


3600 


MOS-3 



"AH the spectra taken during the observation. 

'Out of all the spectra taken, 11 30-min exposures were used to obtain the final combined spectrum. 
°Out of all the spectra taJsen, 22 30-min exposures were used to obtain the final combined spectrum. 
'^We did not use data taken at this night because of their low quality (bad seeing). 
'^Not identified yet. The data analyses are in progress. 
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Table 6 

Spectroscopic properties of the Lyo emission of the z = 6.96 LAE 



Object z 


F(Lya) 
(10-"erg s-l cm-2) 


-L(Lya) 
(lO^erg s-l) 


6'Fi?(Lya) 
(Mo yr-i) 


FWHM 

(A) (kms-l) 


Sw 

(A) 


S/N 


IOK-1 6.96 


2.00 


1.13 


10.24 


13 403 


9.46 ± 0.39 


5.5 



z = l LAE 
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Table 7 

Spectroscopic properties of the UV continuum of the z = 6.96 LAE 



Object z 


(lO^^erg s-i Hz-i) 


SFR{UV} 
(Mo yr-i) 


IOK-1 6.96 


2.58 


36.1 
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Table 8 

The neutral IGM fractions obtained by several independent methods 





Neutral fractions xhi 




Method 


z ~ 6 z 6.3 


z = 6.6 


z = 7.0 


(1) Quasar GP test" 


0.01-0.04 — 






(2) Gamma Ray Burst'' 


— < 0.17-0.60 






(3) Lya LF" 




< 0.45 




(4) Lya LF and LBG UVLF'^ 




- 0.12-0.42 


- 0.12-0.54 


(5) Model and observed Lya LFs-^ 




- 0.24-0.36 


- 0.32-0.64 



fP IFan et all ||20M ). 

(2) Totani ct al. ( 2006). 

(3) ,Kashikawa ct al. (2006b!). 

(4) This study and Yoshida ct al. ( 20061') . iBoiiwens fc lUingworthl (|2006n and lBouwens et al.l (|2006l V 

(5) This study and iKobavashi et al.l I 12'M7|. K07 model). 
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4000 5000 6000 7000 8000 9000 10000 11000 

Wavelength (A) 

Fig. 1. — Filter transmission of the Suprime-Cam broadbands {BVRi'z': thin solid curves) and narrowbands (NB816, NB921: thin solid 
curves and NB973: thick solid curve) used for our photometry. The OH night skylines are also overplotted with thin dashed curve. Our 
NB973 filter is ~ 1.5 times wider in FWHM than other narrowband filters and includes some OH lines in it s bandpass. The mo del spectrum 
energy distribution (SED) of our target, a z = 7 LAE, obtained using a stellar population synthesis model IIBruzual &: Charlotir2003:) with a 
metallicity of Z = Zq = 0.02, an age of f = 1 Gyr, Salpeter initial mass function with lower and upper mass cutoffs of mi, =0.1 Mq and 
mu = 100 Mq and exponentially decaying star formation history for t = 1 Gyr and with a Lya emission of the rest frame equivalent width 
of 50A, is also shown with thick dashed line. The NB973 flux of a 2 = 7 LAE is expected to show a strong excess with respect to its z' band 
flux and should not be detected in other shortward wavebands. 
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0.7 ' 

21.5 22 22.5 23 23.5 24 24.5 25 

NB973 (total) [mag] 

Fig. 2. — The detection completeness of our NB973 image of the SDF, calculated for every 0.5 mag bin. The dashed line shows our survey 
limit. The completeness does not reach 1.0 even for the objects with bright magnitudes. This is because the blended or overlapped objects 
tend to be counted as one object by the SExtractor. The completeness is corrected when the number and luminosity densities of 2 = 7 LAE 
are calculated in §??. 
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Fig. 3.— z'- NB973 vs. i' - z' plot of the objects with NB973 < 24.9 (total mag.) detected in the SDF (shown by black dots). The upper 
right rectangle region surrounded by the solid line indicates the color selection criterion (2). Colors and selection criteria of IOK-1, IOK-2 
(larger green filled squares) and IOK-3 (larger blue filled circle) as well as Obj-4 and Obj-5 (filled diamonds) are also shown and labeled. 
The colors of model LAEs at z = 5.0-7.0 with the rest frame Lyo Une equivalent widths of EWoil^ja) = 0, 10, 20, 50, 100, 150, 200, 250 and 
300A are shown in other several different symbols and lines as labeled in the diagram. Each point with i?Wo(Lya) = 0, which is the first 
point in each sequence, is circled. The model 2 = 7 LAEs are smaller red filled squares with a solid line. All the i' and z' magnitudes fainter 
than their 2a limits were replaced by the 2a values in the application of the color selection criteria and are shown by arrows in this diagram. 
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B V R r NBBIS Z' NBS21NB&73 




Fig. 4. — The muti-waveband images of the IOK-1, IOK-2 and IOK-3 as well as NB973-3(7 excess objects, Obj-4 and Obj-5 (labeled by 
1, 2, 3, 4 and 5, respectively). IOK-1, a. z = 6.96 LAE, and IOK-2 are clearly detected only in NB973. IOK-3, a z = 6.6 LAE identified 
byHCashikawa et al. (2006b), shows a significant excess in both NB921 and NB973 against z' at the same time but is obviously brighter in 
NB921. Obj-4 is seen in all the narrowbands, i' and z' bands while Obj-5 is detected in z' , NB921 and NB973. Both objects show 3(T-excess 
in z'— NB973 but no excess in NB816 (bandpass for z = 5.65—5.75 LAEs), NB921 (bandpass for z = 6.5-6.6 LAEs) and thus could be either 
of 2; = 6.2-6.4 galaxies, low-z EROs or late type dwarf stars. 
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Fig. 5. — z'— NB973 (2" aperture mags.) color as a function of NB973 (total) magnitude of the objects detected in the SDF (shown by 
dots). The dotted curve shows 3cr error track of z'— NB973 color. The horizontal solid line is a part of our color selection criterion (2), z' — 
NB973 > 1.0. The vertical dashed line indicates the detection limiting magnitude of our survey, NB973 = 24.9(5(t). The diagonal dashed line 
is the 2cr limits of i' and z' 2" aperture magnitudes. The IOK-1, -2 and -3 as well as Obj-4 and -5 are denoted by the same symbols as those 
in Figure |3] 
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N(i') - X 
NB973(5a) = 24.9 
i'(5c7) - 0.33 = 26.52 




1000 

22.5 23 23.5 24 24.5 25 25.5 26 26.5 27 

NB973 = i' - 0.33 [total mag] 

Fig. 6. — Comparison of the number counts of objects detected in the NB973 (plus symbols) and i' band (crosses) images of SDF down 
to their detection limits, NB973 = 24.9 (Sir) (dashed line) and i' = 26.85 (5it) (dotted line), respectively. We try to extrapolate the NB973 
number count of objects between NB973 = 24.9 and 26.52 by using the mean color relation of < i'— NB973>= 0.33 and i' band number 
count at the corresponding magnitude range. Some fraction of such objects can be candidate variables that would become brighter than 
NB973 = 24.9 at another occasion. Since i' band number count is slightly (xl.l-xl.2) higher than that of NB973 at NB973 < 24.9, our 
extrapolation can be an overestimation. 
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Fig. 7. — Cumulative hya LFs of LAEs at z = 5.7 (Shimasaku e t ani2006l . filled squares for their photmetric LF), 6.6 l IKashikawa et aTl 
I2006bl . filled circles and triangles for their photometric (upper limit) and spectroscopic (lower limit) LFs, respectively) and 7 (this study, the 
diamond for IOK-1 spectroscopic data). All the errors include cosmic variance and Poissonian errors for small-number statistics (See US. II 
for details). All the data and errors are corrected for their detection completeness. The long-dashed, short-dashed and solid curves are the 
intrinsic (i.e., not affected by neutral IGM) Ly« LFs at z = 5.7, 6.6 and 7, respectively, predicted by K07 LAE evolution model (See i|5.2.3ll . 
The vertical dashed line shows our survey limit to which the LFs are integrated down to obtain LAE number and luminosity densities. 
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Age of the Universe (Gyr) 

0.97 0.810.75 




5 

Redshift 

Fig. 8. — The number density ni^ya, Lyo line luminosity density pLya a-^d star formation rate density SFRDi^yc of LAEs at z = 5.7, 6.6 and 
7 derived from the latest Subaru/Suprime-Cam LAE surveys and those at 2.3 < z < 5.7 from liter ature down to Lii^ii^fLy o:) = 1.0 X 10^^ erg 
s~^. The densities at z = 5.7, 6.6 and 7 i n the SDF are calc ulated from the photometric sample of lShimasaku et al.l I I2006I . open squares), the 
photometric and spectroscopic samples of lKashikawa et al.l |j2006b, open circles and triangles, respectively), and IOK-1 s pectrum (large filled 
circles ), respectively. Densities a t 2.3 < z < 4.5 and z ^ 4.5 are calculated using the best- fit Lya Schechter LPs from Ivan Breukelen et aU 
l|2005l . inverse triangle) and Daws on et al.l 1)20071 . pentagon), respectively . Also , densities at z = 3.1, 3.7 and 5.7 in the ~ 1.0 deg^ of SXDS 
field are calculated using the best-fit Lya Schechter LPs from lOuchi et ahl II2007I . diamonds) Each horizontal error bar shows the redshift range 
of each survey. The vertical error bars at z = 5.7, 6.6 and 7 include both cosmic variance and Poissonian errors for small-number statistics 
while those at z < 5.7 contain only cosmic variance since Poissonian errors for them are negligibly small (See i|5.1l for details). The data and 
vertical error bars at z = 5.7, 6.6 and 7 are corrected for their detection completeness. The plus symbols at z = 5.7, 6.6 and 7 with solid 
lines show the expected densities obtained by integrating the intrinsic (i.e, not affected by neutral IGM) Lyo LPs predicted by the K07 LAE 
evolution model in i|5.2.3l At z > 5.7, the densities clearly decrease with increasing redshifts and smaller than the model-predicted values, 
implying that the Lya lines might be attenuated by the possibly increasing neutral IGM at the reionization epoch. 
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z=5.7 Hu et al. 2004 
z=5.7 Shimasaku et al. 2006 
z=6.6 Kashikawa et al. 2006 
z=7.0 This study 
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Fig. 9. — The rest frame UVLFs (per unit absolute magnitude) of 2 = 5.7 and 6.6 LAEs compared with that at 2 = 7 derived using IOK-1 
data. Our survey detection hmits at 5cr and 3cr are shown by vertical dashed lines. The errors include cosmic variance and Poissonian errors 
for small-number statistics (See i|5.1l for details). All the data and errors are corrected for their detection completeness. The 2 = 7 UVLF 
does not evolve from 2 = 5.7-6.6. The long-dashed, short-dashed and solid curves are UVLFs of LAEs at 2 = 5.7, 6.6 and 7 predicted by K07 
model, respectively. They are in good agreement with the observation data. 
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Fig. 10. — Cumulative Lyo LFs of LAEs at 2 = 5.7 (Shimasaku et al.l'20i)(3, filled squares) and 7 (this study; diamond). Errors include 
cosmic variance and Poissonian errors for small- number statistics (Sec ^S.ll for details). All the data and errors are corrected for their detection 
completeness. The solid curve is the best-fit z = 5.7 Schechter LF (See ^5.11 1. The dashed curve is the same Lya Schechter LF but made to 
evolve fr om z = 5.7 to 2 = 7 b y imposing /^*><=^P<='^* ^ ^*=5 7 ^ 0.58 obtained using the correlation of AM^y / Az ~ 0.47 extrapolated from the 
result of lYoshida et all II2006I . See their Figure 22) and assuming the same correlation can also hold for Lyo luminosity (See text in i|5.2.2l l. 



